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Diffraction of light sets the fundamental limit for optical lithography. Many quantum lithography 
schemes have so far been proposed to overcome this limit either by making use of highly entangled 
photons, multi-photon processes or multiple Lambda transitions, which are all experimentally high- 
demanding. Recently, Liao et al. proposed a novel quantum lithography scheme which merely 
employs Rabi oscillation to surpass the diffraction limit. Here we report a faithful experimental 
realization of this scheme. Resolution up to ninth of the Rayleigh diffraction limit has been observed. 

Possibility of creating an arbitrary pattern is also tested experimentally by demonstrating the peak 
narrowing process using several Rabi hoppings together with state-selective optical depletion. Our 
work may have direct applications in atom pattern engineering for quantum information or quantum 
simulation applications, and will also possibly boost the adoption of quantum lithography into real- 
world applications in the near future. 

PACS numbers: 42.50.St, 42.50.Gy, 42.50.Ct, 85.40.Hp 


Optical lithography is the most important method to 
fabricate microscale structures in modern semiconductor 
industry jl]. The resolution of optical lithography is fun¬ 
damentally limited by optical diffraction pQ. Many quan¬ 
tum lithography schemes [2] have been proposed so far to 
create patterns with both the feature size and feature sep¬ 
aration overcoming the diffraction limit, by either making 
use of the quantum nature of light or the quantum na¬ 
ture of matter. By making use of highly path-entangled 
NOON states 3j [4], one can create interference patterns 
with a resolution which is enhanced by N times over the 
diffraction limit. Nevertheless, high-order NOON states 
with high intensities are extremely difficult to prepare. It 
is also possible to create sub-wavelength features by mak¬ 
ing use of multi-photon nonlinear processes with classical 
light [5- HE], however the nonlinearity is typically rather 
small and extremely high laser intensities are required. 
There is also proposal which makes use of dark states in 
a multiple Lambda system and requires only moderate 
laser intensity [9], nevertheless this scheme suffers from 
the difficulty of finding the well-structured multi-Lambda 
energy system and the adoption of many laser beams. 

Recently, Liao et al. have proposed a novel quan¬ 
tum lithography scheme m, which simply involves two 
atomic levels and can in principle reduce the feature 
size and separation beyond the diffraction limit arbitrar¬ 
ily. This scheme relies on the process of Rabi oscillation 
driven by a standing-wave light field. The standing-wave 
field can be created by two-beam interference, thus giv¬ 
ing an intensity period of A e ff = A/2sin 0, where 0 is 
half of the separation angle between the beams. As the 
Rabi frequency depends on the electric field strength, 
UIr = d • E/h, nonlinearity of the light field is directly 


transferred to the spatial distribution of CLr. At posi¬ 
tions with maximal electric field strength, atoms oscil¬ 
late fastest between internal states; while at positions of 
interfering minima, which ideally is zero, atoms can not 
be excited. If a n pulse (for the peak UIr) is applied, 
an atomic pattern with periodic spacing of A e ff can be 
created. If a Ntt pulse is applied, the atomic pattern 
resolution can be reduced to X e ^/N, which in principle 
has no upper limit. As a contrast, in traditional opti¬ 
cal imaging and optical lithography, the Rayleigh crite¬ 
rion puts an upper limit of about X/2N.A. in resolution, 
where N.A. is the numerical aperture. For lithography 
in air, this limit is A/2 sin 0, where 0 here is half of the 
included angle of the objective lens. Thus, with the Rabi 
oscillation method, by applying a Ntt pulse, resolution 
can be enhanced by N times over the diffraction limit in 
classical optical lithography. 

In this paper, we report the first experimental real¬ 
ization of this scheme. Two-photon Raman transitions 
are employed to induce Rabi oscillations between two in¬ 
ternal ground levels. One Raman beam is periodically 
patterned through interference, which gives A e ff ~ 57/im, 
while the other is flat for simplicity. We manipulate laser- 
cooled atoms with these Raman beams, and make use of 
a specially designed absorption imaging system to mea¬ 
sure the atomic patterns. High-resolution spatial pat¬ 
terns up to ninth of the Rayleigh diffraction limit have 
been created and observed for atoms in specific hyper- 
fine state. Moreover, by combining the Rabi oscillation 
together with state-selective optical depletion, arbitrary 
patterns can be produced in principle m- In our ex¬ 
periment, we demonstrate the process of peak narrowing 
without changing the peak separation by applying mul- 
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tiple pulses. 




FIG. 1: Experimental scheme, (a) Energy levels. The D 1 
line transition of 87 Rb is used as the channel for two-photon 
stimulated Raman transitions, as well as for i r pumping. But 
for absorption imaging, F — 2 —»> F' — 3 cycling transition 
of D2 line is used, (b) Illustration of spatial distribution of 
Rabi frequency and atomic pattern on each hyperhne ground 
state, (c) Typical atomic patterns in experiment. Left: in the 
dipole trap; Middle: after a In Raman pulse. Right: after a 
2n pulse, (d) Experimental layout. A cold atomic ensemble 
is captured with a dipole trap inside the vacuum glass cell. 
Then a simple interferometer with two arms is placed outside 
the vacuum, where Raman beams are splited and polarization 
rotated, then interfere at the center. After the cell, a high 
resolution lens group compensates the abberation introduced 
by the glass cell and then a 780 nm narrow bandpass filter 
(B.F.) is placed before the CCD camera to filter out other 
scattered photons except the imaging beam. 

The detailed experimental scheme is shown in Fig.[^ A 
cold 87 Rb atomic ensemble is prepared by first capturing 
with a standard magneto-optical trap (MOT) and then 
a molasses cooling process which results in an atom tem¬ 
perature of - 10 fi K. Afterwards, the atoms are loaded 
into an optical dipole trap m to compress the radial 
dimension of the atomic ensemble. The dipole trap is 
provided by a tightly focused 850 nm laser beam from a 
diode laser, with a waist of about 30 fi m and 170 mW 
power, thus giving a trapping depth of about 57 fiK. 
The atoms held in the dipole trap have a width of about 
27 /im and peak optical depth of about 0.3, as shown in 
Fig. B c). Before driving Rabi oscillations, these atoms 
are released from the dipole trap, and then optically 
pumped to the state of \F = l,rriF = 0). To image 
the projected atom patterns down to several fi m, we use 
a high resolution imaging design similar to Ref. [12], to 
compensate the spherical aberrations from the vacuum 
glass cell. The imaging system has a magnification of 9.7 
and a resolution better than 4.5 yum, as measured with 
a Thorlabs test target. An Andor IXon EMCCD, with 
pixel size of 16 fi m, is used to record images for off-line 
analysis. Absorption imaging is taken with the D2 line 
cycling transition of 87 Rb atoms, where only atoms on 


F = 2 ground state are detected. 

To create the standing-wave light fields for driving 
Rabi oscillations, we modify the original two-level scheme 
to a three-level stimulated Raman transition process m- 
As shown in Fig. [lja), two hyperfine ground sublevels 
are connected via a coherent two-photon transition pro¬ 
cess, where the Raman fields are provided by two phase 
locked diode lasers, and each with a single photon de¬ 
tuning much larger than the nature linewidth of D1 line, 
A ^ T. For 87 Rb atoms, there are two possible paths 
for two-photon transitions between | F = 1 ,ijif = 0) and 
|F = 2, mi? = 0) ground sublevels. Given two Raman 
fields with H and V polarizations, constructive interfer¬ 
ence will happen between these two paths due to the 
Clebsch-Gordan coefficients of relative transitions, thus 
two-photon transitions can still be driven between these 
two sublevel states. 

The layout of the Raman fields is shown in Fig. 0 d )> 
where one of the Raman beams is splitted into two paths 
by a Glan-Taylor prism, and then the polarization of 
the reflected path is rotated back to be horizontal to 
make them interfere at the center, leaving a standing- 
wave light field there. The other beam, however, only 
transmits through one path of the prism and forms an 
homogeneous field at the center. The two-photon Ra¬ 
man Rabi frequency depends on single-beam Rabi fre¬ 
quencies, Qr(x) oc ck Ei(x)E 2 (x ), where 

fli(x) = d • Ei(x)/h and i = 1,2 correspond to light 
fields of Raman beam 1 and 2. Note the ideal inter¬ 
fering field is Ei(x) = 2E® cos7nr/A e ff and the flat field 
is E 2 (x) = E ®, thus Raman Rabi frequency can be de¬ 
scribed as CIr(x) = Q° r cos7nr/A e ff, where A e ff ~ 57 yum 
for half separation angle of 0 = 0.4°. Note if both beams 
are splitted into two paths and then interfere with them¬ 
selves at the center, relative path length control of the 
interfering arms is necessary to make sure both interfer¬ 
ing patterns are in-phase or with a special phase offset. 
And due to the simplified configuration here, atoms on 
F = 1 state experience position dependent light shifts, 
A 1 ac { x ) c* while atoms on F = 2 have homoge¬ 

neous shifts. This spatial periodic differential light shifts, 
created along the axial direction of atoms in the dipole 
trap, will limit the spatial visibility at some positions but 
will not affect the resolution enhancement according to 
our simulations. 

In our experiment, the Raman beams have a diam¬ 
eter of 1.3 mm, each has a single photon detuning of 
A = +750 MHz relative to F f = 2 state of the D1 line. 
Raman beam 1 has a total power of 1 mW, while beam 
2 has a power of 2 mW to match the light shifts for 
atoms at interference peak positions. With these beams, 
the on resonance Raman Rabi frequency at peak posi¬ 
tions is calculated to be about 2n x 550 kHz. The two- 
photon detuning frequency is set to be 240 kHz, which 
is to match the on resonance condition for atoms at half 
of the maximum Rabi frequency positions. Taking ac- 
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FIG. 2: Resolution enhancement with various Raman pulse 
durations in a single A e ff range. (a),(b) Simulated and exper¬ 
imental atomic patterns in a single A e fr range. The measured 
pictures are interpolated with 4 points in each pixel area to 
show the patterns clearly. It’s demonstrated with a pulse 
duration of 97r, up to 9 peaks can be created in the single 
intensity period area, (c) Fitted pattern visibility for several 
pulse durations, where the horizontal axis is position in space 
with unit of /im, and vertical axis is optical depth. These 
points used in fitting visibility are raw pixel points. Visibili¬ 
ties for l7r, 27r, 47r and 57r durations, are 0.83, 0.79, 0.48, and 
0.39 as fitted, respectively. 


mainly contributed by the relative phase fluctuations be¬ 
tween two interfering arms. Second, atoms are heated 
up during releasing from the trap and optical pumping. 
Diffusion of atoms will finally erase the atomic patterns 
for several 10 fi s, thus the patterns have to be imaged 
immediately after creation. 
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count of detuning, the spatial maximum Rabi frequency 
becomes 2tt x 605 kHz, which agrees well with our mea¬ 
surements. In Fig. [ljc), periodic atomic patterns clearly 
appear after applying l7r or 2tt pulses with duration of 
0.8 fi s or 1.6 /is, respectively. Considering the pixel size 
and magnification ratio, the periodic separation is fitted 
to be Ag ff = 59.3 fi m for l7r pulse, which corresponds to 
the diffraction limit if projecting an image through this 
small aperture. 

After the initial exploration with short tt pulses in 
Fig. B c), we now turn to multi- tt pulses. To present 
the structures even clearer, we now focus on the area of 
a single effective wavelength range, as shown in Fig. [2|a- 
b). It clearly demonstrates that up to nine peaks can be 
created in a single A e ff range, where the smallest sepa¬ 
ration reaches down to about 6.5 fi m. For short pulse 
durations up to 57r, we find it’s able to fit the pattern 
distribution with normal sinusoidal functions very well 
as shown in Fig. [2jc) , and the visibility is calculated as 
the amplitude of oscillation divided by the average. For 
longer pulse durations, although it’s able to identify num¬ 
ber of peaks, the visibility gets reduced severely. This is 


FIG. 3: Narrowing of the atomic pattern with different pulse 
combinations, (a) Sequence of the experiment. Atomic pat¬ 
terns are first created on F = 1 state, and then narrowed with 
the second Raman pulse. Then after all atoms originally on 
F — 2 state are pushed away of the area, the patterns are 
repumped to F — 2 for imaging, (b) Experimental results of 
width narrowing, with the first pattern being created by a 7r 
pulse. Left, absorption imaging of the cloud, images here are 
rotated relative to in Fig. [2] Right, cross-section view of the 
axial width, where horizontal axis denotes the position in real 
space with unit of /im and vertical axis denotes optical depth. 
For specific pulse combinations and peak widths, see in text. 

With only a single Raman pulse, we have already 
shown the diffraction limit can be easily circumvented via 
Rabi oscillations driven by a standing-wave field. With 
a combination of several pulses, there are more excit¬ 
ing prospects with this method, i.e., creating arbitrary 
patterns as proposed in the original scheme m- Here 
we show how to create special patterns with two pulses, 
i.e. narrowing the peak width but without altering the 
separation distance, as shown in Fig. [3j All atoms are 
initially prepared on |F = l^mp = 0) state, then a tt 
pulse is applied to create the initial atomic pattern on 
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F = 1 state. Subsequently, a strong depleting beam res¬ 
onant with the D 2 line cycling transition is applied in 
vertical direction during 20 fi s, thus pushing all atoms 
on F = 2 state out of the imaging area. Afterwards, the 
second Raman pulse with specific duration is applied, to 
narrow the peak width of atoms on F = 1 state. Then, 
the pushing beam with 20 fi s duration is shined again to 
clear all atoms transferred to F = 2 state. In the end, a 
repumping beam is applied to transfer the spatial distri¬ 
bution from F = 1 to F = 2 state for absorption imaging. 
The experimental results are shown in Fig. [3^b), where 
the second Rabi pulse with durations of zero, l7r, |7r, 
27r, |7r are used. For each duration, 20 pictures are av¬ 
eraged for analysis. The fitted widths for these pulses 
are found to be 38.4 /i m, 37.8 /im, 31.2 /i m, 22.3 /im, 
17.8 /im, and 14.5 /im, where the smallest width realized 
here is about 1/4 of the diffraction limit. Note that, be¬ 
fore taking these images here, atoms on F = 1 state 
have already thermally diffused, during pushing away 
and repumping processes. As a reference, the theoret¬ 
ical predicted widths for these pulse combinations are, 
31.0 /im, 24.6 /im, 18.5 /im, 14.0 /im, and 12.8 /im, re¬ 
spectively, where the slightly larger widths measured can 
be explained as above. 

There are several technical imperfections yet to be re¬ 
solved in future works. First, with traditional optical 
imaging, the detection itself is diffraction limited, that’s 
why we have to explore the diffraction nature at scales 
much larger than optical wavelength. To fully explore the 
ultimate performance of the scheme for grazing incident 
Raman beams, at the scale of subwavelength, other meth¬ 
ods like atomic force microscope have to be utilized to 
measure after imprinting the pattern on a substrate. Sec¬ 
ond, relative phase of the interference arms hasn’t been 
stabilized in this work. Without phase control, the inter¬ 
ference field randomly drifts at the center at the scale of 
minutes, which is the same time scale of taking pictures 
for averaging in the experiment. This limits the maximal 
resolution can be achieved in our current work. Third, to 
simplify analysis and control, only single Raman beam is 
splitted and then interfering at the center. This results 
in inhomogeneous differential Stark shifts between inter¬ 
nal levels across the standing-wave area, and limits the 
uniformity of visibility across the sample. The inhomo¬ 
geneous shifts can be cancelled by making both of the 
Raman beams interfere at the center as well as making 
sure both standing waves are in-phase with each other. 
Finally, to create arbitrary patterns as proposed in the 
original scheme, multiple Raman pulses with different du¬ 
rations, controllable node positions of the standing wave, 
as well as a continuous incoming atomic source are nec¬ 
essary. 

To summarize, the experiment presented here provides 
a clear demonstration of the quantum lithography scheme 
based on Rabi oscillations, and shows that resolution 
can be enhanced by up to nine times over the diffrac¬ 


tion limit. With multiple Raman pulses, we can also 
create much narrower atomic peaks without altering the 
peak separations, where the smallest width is about one- 
fourth of the diffraction limit. Although patterns cre¬ 
ated are only one dimensional here, it’s straightforward 
to extend this scheme to two-dimension, only by using 
the other Raman beam to create a standing-wave field 
along the other direction, to create a two-dimensional 
Rabi frequency field. In comparison with the traditional 
atom lithography experiments IMS! or recent atom lo¬ 
calization experiments IMS] with dark states, our ex¬ 
periment highlights the advantage that both the feature 
size and feature separation can overcome the diffraction 
limit. Our work may have applications in optical lithog¬ 
raphy, atom pattern engineering for quantum informa¬ 
tion or quantum simulation experiments, and could also 
possibly boost the adoption of quantum lithography into 
real-world applications in the near future. 

We thank M. Suhail Zubairy for stimulating dis¬ 
cussions. This work was supported by the National 
Natural Science Foundation of China, National Funda¬ 
mental Research Program of China (under Grant No. 
2011CB921300), and the Chinese Academy of Sciences. 
X.-H. B. and B. Z. acknowledge support from the Youth 
Qianren Program. 


[1] T. Ito and S. Okazaki, Nature 406, 1027 (2000). 

[2] M. Al-Amri, Z. Liao, and M. S. Zubairy, in Advances in 
Atomic, Molecular, and Optical Physics , edited by E. A. 
Paul Berman and C. Lin (Academic Press, 2012), vol. 61 
of Advances In Atomic, Molecular, and Optical Physics , 
pp. 409 - 466. 

[3] A. N. Boto, P. Kok, D. S. Abrams, S. L. Braunstein, C. P. 
Williams, and J. P. Dowling, Phys. Rev. Lett. 85, 2733 
( 2000 ). 

[4] M. D’Angelo, M. Chekhova, and Y. Shih, Phys. Rev. 
Lett. 87, 013602 (2001). 

[5] G. S. Agarwal, R. W. Boyd, E. M. Nagasako, and S. J. 
Bentley, Phys. Rev. Lett. 86, 1389 (2001). 

[6] P. Hemmer, A. Muthukrishnan, M. O. Scully, and 
M. Zubairy, Phys. Rev. Lett. 96, 163603 (2006). 

[7] S. J. Bentley and R. W. Boyd, Opt. Express 12, 5735 
(2004). 

[8] P. Avi, B. Dayan, M. Vucelja, Y. Silberberg, and A. A. 
Friesem, Opt. Express 12, 6600 (2004). 

[9] M. Kiffner, J. Evers, and M. S. Zubairy, Phys. Rev. Lett. 
100, 073602 (2008). 

[10] Z. Liao, M. Al-Amri, and M. Suhail Zubairy, Phys. Rev. 
Lett. 105, 183601 (2010). 

[11] S. J. M. Kuppens, K. L. Corwin, K. W. Miller, T. E. 
Chupp, and C. E. Wieman, Phys. Rev. A 62, 013406 
( 2000 ). 

[12] J. Mun, Ph.D. thesis, MIT, Boston (2008). 

[13] D. S. weiss, M. A. Kasevich, and S. Chu, Phys. Rev. A 
45, 342 (1992). 

[14] G. Timp, R. E. Behringer, D. M. Tennant, and J. E. 
Cunningham, Phys. Rev. Lett. 69, 1636 (1992). 



5 


[15] K. S. Johnson, J. H. Thywissen, N. H. Dekker, K. K. 
Berggren, A. P. Chu, R. Younkin, and M. Prentiss, Sci¬ 
ence 280, 1583 (1998). 

[16] M. Miitzel, S. Tandler, D. Haubrich, D. Meschede, 
K. Peithmann, M. Flaspohler, and K. Buse, Phys. Rev. 
Lett. 88, 083601 (2002). 

[17] H. Li, V. Sautenkov, M. Kash, A. Sokolov, G. Welch, 


Y. Rostovtsev, M. Zubairy, and M. O. Scully, Phys. Rev. 
A 78, 013803 (2008). 

[18] N. A. Proite, Z. J. Simmons, and D. D. Yavuz, Phys. 
Rev. A 83, 041803 (2011). 

[19] J. A. Miles, Z. J. Simmons, and D. D. Yavuz, Phys. Rev. 
X 3, 031014 (2013). 



